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Abstract 

This  thesis  examines  the  effects  of  multipath  interference  on  Low  Probability 
of  Intercept  (LPI)  waveforms  generated  using  intersymbol  dither.  LPI  waveforms 
are  designed  to  be  difficult  for  non-cooperative  receivers  to  detect  and  manipulate, 
and  have  many  uses  in  secure  communications  applications.  In  prior  research,  such 
a  waveform  was  designed  using  a  dither  algorithm  to  vary  the  time  between  the 
transmission  of  data  symbols  in  a  communication  system.  This  work  showed  that  such 
a  method  can  be  used  to  frustrate  attempts  to  use  non-cooperative  receiver  algorithms 
to  recover  the  data.  This  thesis  expands  on  prior  work  by  examining  the  effects  of 
multipath  interference  on  cooperative  and  non-cooperative  receiver  performance  to 
assess  the  above  method’s  effectiveness  using  a  more  realistic  model  of  the  physical 
transmission  channel. 

Both  two  and  four  ray  multipath  interference  channel  models  were  randomly 
generated  using  typical  multipath  power  profiles  found  in  existing  literature.  Differ¬ 
ent  combinations  of  maximum  allowable  symbol  delay,  pulse  shapes  and  multipath 
channels  were  used  to  examine  the  bit  error  rate  performance  of  1)  a  Minimum  Mean 
Squared  Error  (MMSE)  cooperative  equalizer  structure  with  prior  knowledge  of  the 
dither  pattern  and  2)  a  Constant  Modulus  Algorithm  (CMA)  non-cooperative  equal¬ 
izer.  Cooperative  MMSE  equalization  resulted  in  approximately  6-8  dB  BER  perfor¬ 
mance  improvement  in  Eh/N0  over  non-cooperative  equalization,  and  for  a  full  range 
symbol  timing  dither  non-cooperative  equalization  yields  a  theoretical  BER  limit  of 
T(,=10-1.  For  50  randomly  generated  multipath  channels,  six  of  the  four  ray  chan¬ 
nels  and  15  of  the  two  ray  channels  exhibited  extremely  poor  equalization  results, 
indicating  a  level  of  algorithm  sensitivity  to  multipath  conditions. 
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Low  Probability  of  Intercept  Waveforms 
via  Intersymbol  Dither 
Performance  Under  Multipath  Conditions 

I.  Introduction 

1 . 1  Motivation 

As  the  role  of  technology  grows  on  the  battlefield,  the  ability  to  sustain  se¬ 
cure  communications  links  has  become  vital  to  military  operations  [9].  With 
this  need  for  secure  communications  over  a  broad  range  of  activities  and  environ¬ 
ments,  wireless  communications  systems  have  become  increasingly  popular  for  mili¬ 
tary  use  [15].  The  flexibility  provided  by  wireless  broadcasting  systems  transmitting 
over  large  areas  also  makes  them  vulnerable  to  detection  and  eavesdropping  by  third 
parties.  As  such,  a  great  deal  of  current  communications  research  focuses  on  ways  to 
mitigate  the  risk  of  eavesdropping,  detection,  and  exploitation  of  wireless  communi¬ 
cation  systems. 

Three  basic  areas  of  research  into  reducing  the  vulnerability  of  friendly  trans¬ 
missions  to  eavesdropping  and  exploitation  are  Low  Probability  of  Detection  (LPD), 
Low  Probability  of  Interception  (LPI),  and  Low  Probability  of  Exploitation  (LPE). 
LPD  and  LPI  deal  with  methods  to  prevent  or  make  it  difficult  for  a  non-cooperative 
receiver  to  be  able  to  pick  out  the  transmitted  signal  of  interest  from  environmental 
noise  or  other  signals  [14].  For  example,  spread  spectrum  and  frequency  hopping  are 
LPI/LPD  methods  in  that  they  make  the  transmitted  signal  difficult  to  detect  or 
intercept  [17].  LPE  methods  focus  on  preventing  any  meaningful  extraction  of  data 
from  the  transmitted  signal  of  interest  by  non-cooperative  receivers  and  third  parties. 

Multipath  interference  is  a  common  issue  in  wireless  communication  systems, 
and  is  caused  by  a  transmitted  signal  reaching  the  receiver  by  one  or  more  time 
delayed  paths  by  reflecting  off  of  objects  in  the  transmission  environment  [14].  These 
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delayed  transmission  paths  can  then  interfere  with  the  desired  transmitted  signal, 
causing  a  degredation  in  system  performance.  This  thesis  examines  the  effects  of 
multipath  interference  on  the  intentional  Intersymbol  Interference  (ISI)  waveform 
designed  in  [8]  to  prevent  non-cooperative  receivers  from  determining  the  transmitted 
communications  symbols. 

1.2  Problem  Statement 

Because  of  the  proliferation  and  widespread  use  of  wireless  communications  for 
military  and  civilian  use,  communications  have  become  vulnerable  to  eavesdropping, 
interception,  and  exploitation  by  unintended  recipients.  As  such,  strategies  intended 
to  reduce  these  risks  associated  with  wireless  broadcast  transmissions  such  as  LPI, 
LPD,  and  LPE  techniques  are  desirable.  One  strategy  is  to  utilize  techniques  designed 
to  make  it  difficult  for  non-cooperative  receiver  designs  to  detect  or  intercept  the 
communications  waveform.  Previous  research  by  [8]  developed  a  waveform  designed 
to  perform  this  task,  and  validated  its  effectiveness  over  a  Line  of  Sight  (LOS)  channel. 

This  thesis  will  develop  several  different  multipath  models,  validate  them,  and 
test  the  performance  of  the  ISI  waveform  developed  in  [8]  under  multipath  conditions. 
The  ISI  waveform  modifies  Differential  Phase  Shift  Keyed  (DPSK)  modulation  to  frus¬ 
trate  non  cooperative  receiver  strategies.  This  is  accomplished  by  inserting  ISI  into 
the  transmitted  waveform  to  disrupt  accurate  non  cooperative  symbol  estimation  [8] . 
The  ISI  waveform  achieved  symbol  estimation  accuracy  close  to  optimal  DPSK  per¬ 
formance  over  a  LOS  channel  at  the  cost  of  additional  receiver  complexity  to  mitigate 
the  ISI.  Any  additional  complexity  incurred  by  the  addition  of  multipath  interference 
is  examined,  as  well  as  the  potential  performance  loss  of  a  non-cooperative  receiver 
compared  to  a  cooperative  receiver  designed  to  account  for  ISI  and  multipath. 

1 . 3  Thesis  Organization 

Chapter  2  reviews  the  pertinent  background  information  to  the  problem,  and 
provides  any  necessary  conceptual  derivations.  Topics  covered  include  multipath 
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Interference,  the  design  of  the  ISI  waveform,  and  the  communications  model  used 
throughout  the  thesis.  Chapter  3  details  the  development  of  the  models  and  simu¬ 
lations,  and  provides  validation  of  these  models.  This  includes  the  selection  of  the 
multipath  models,  and  the  development  and  design  of  a  cooperative  receiver  model 
to  equalize  the  ISI  and  multipath.  Chapter  4  presents  the  simulation  results,  detail¬ 
ing  the  performance  of  the  waveform  under  multipath  conditions.  Comparisons  to 
analytical  results  and  previous  research  are  also  presented.  Chapter  5  presents  the 
conclusions  drawn  from  simulations  in  Chapter  4,  and  provides  recommendations  for 
future  research  in  the  problem  area. 
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II.  Background  and  Literature  Review 

2.1  Overview 

This  section  begins  by  providing  an  overview  of  some  common  strategies  for 
securing  wireless  communications.  The  remainder  of  the  chapter  covers  back¬ 
ground  related  to  the  modeling  of  ISI  in  wireless  communications.  A  review  of  types 
of  multipath  models  is  presented,  along  with  the  different  aspects  of  multipath.  Two 
basic  classes  of  equalizers  are  presented;  cooperative  and  non-cooperative,  or  blind 
equalizers.  Cooperative  equalizers  make  use  of  some  apriori,  or  prior,  knowledge  of 
the  transmitted  signal  of  interest  in  order  to  equalize  ISI  and  other  degrading  effects, 
while  blind  equalizers  attempt  to  equalize  these  same  effects  with  minimal  or  no  prior 
knowledge  of  the  transmitted  signal.  An  overview  of  the  intentional  ISI  waveform 
using  symbol  dither  developed  in  [8]  is  also  provided. 

2.2  Existing  Secure  Communications  Strategies 

Encryption  is  one  of  the  most  widely  used  methods  of  securing  communica¬ 
tions,  and  is  used  in  applications  such  as  computer  networking,  secure  mobile  phones, 
etc.  The  main  goal  of  encryption  is  to  prevent  an  unauthorized  user  from  extracting 
meaningful  information  from  the  transmitted  signal,  protecting  the  transmitted  data 
contained  within.  Encryption  techniques  typically  use  an  encryption/decryption  key 
that  is  provided  to  the  intended  users,  and  without  this  decryption  key  a  cooperative 
receiver  will  not  be  able  to  extract  the  data  from  a  transmitted  signal.  The  security 
of  this  method  is  derived  from  the  inability  of  an  unauthorized  party  to  obtain  the 
key,  even  if  they  have  some  knowledge  of  the  encryption  and  transmission  methods 
used  [17]. 

Frequency  hopping  involves  changing  the  carrier  frequency  of  the  modulated 
signal,  typically  according  to  a  pseudorandom  sequence.  A  cooperative  receiver  uses 
knowledge  of  this  pseudorandom  sequence  to  know  when  to  switch  receive  frequencies 
to  locate  and  demodulate  the  signal.  Two  types  of  frequency  hopping  systems  are 
slow  and  fast  hopping  systems.  Slow  hopping  systems  remain  at  a  particular  carrier 
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frequency  for  a  block  of  transmitted  symbols,  while  a  fast  hopping  system  changes 
carrier  frequencies  after  every  symbol.  Without  access  to  the  pseudorandom  sequence 
driving  the  carrier  frequency  changes,  a  non-cooperative  receiver  is  forced  to  locate 
the  new  carrier  after  each  hop  before  attempting  demodulation.  As  such,  a  frequency 
hopping  system  must  switch  carriers  quickly  enough  to  prevent  a  non-cooperative 
receiver  from  being  able  to  detect  the  change  and  match  frequencies  [17]. 

Spread  spectrum  techniques  provide  security  by  making  it  difficult  for  an  unin¬ 
tended  user  to  detect  the  signal  at  all.  Spread  spectrum  methods  smear  the  modulated 
signal  over  a  very  large  bandwidth  by  mixing  it  with  a  pseudorandom  code.  As  this 
disperses  the  signal  power  across  this  large  bandwidth,  it  becomes  much  less  likely 
that  a  third  party  receiver  would  be  able  to  detect  the  transmitted  signal  within  the 
noise  of  the  transmission  environment  [17].  As  in  frequency  hopping  systems,  a  coop¬ 
erative  receiver  must  have  knowledge  of  the  pseudorandom  mixing  code  to  despread 
the  transmitted  waveform  and  reconstruct  the  original  transmitted  signal  for  demod¬ 
ulation.  As  a  non-cooperative  receiver  is  unlikely  to  have  access  to  the  code  used 
for  spreading  the  transmitted  signal  it  is  less  likely  to  be  able  to  extract  any  useful 
information,  providing  an  additional  layer  of  communication  security  [1]. 

2.3  IS  I  and  Equalization 

Digital  communications  systems  transmit  information  as  a  series  of  communi¬ 
cations  symbols  representing  binary  data.  In  a  very  basic  signalling  scheme,  these 
communications  symbols  are  transmitted  back  to  back,  with  the  end  of  one  trans¬ 
mitted  symbol  marking  the  beginning  of  the  next.  As  all  practical  pulse  shapes  used 
for  transmission  of  symbols  are  bandlimited,  there  is  a  limit  to  the  speed  at  which 
the  communications  system  can  alter  the  transmitted  waveform.  This  can  cause  the 
symbols  to  mix  or  blur  together,  creating  ISI  [5] .  Also,  as  these  symbols  are  transmit¬ 
ted  through  some  kind  of  physical  media,  e.g.  through  the  air  in  the  case  of  wireless 
transmissions,  the  received  signal  is  an  attenuated  and  distorted  version  of  the  original 
transmitted  symbol.  These  effects  are  worsened  in  wireless  communications  as  each 


5 


x(n) 


Noise 


x^ct\t) 


b^ct\t) 

-4 


Tsym/O 


y{ct\t) 


Channel 


Sample 


y'[k) 


y\k) 


w' 

1° 

x'(k) 

Filter 

Downsample 

Figure  2.1:  Baseband  communications  model  [8] 


transmitted  symbol  can  take  multiple  paths  to  the  receiver  with  each  path  causing 
independent  amplitude  and  phase  distortions  on  the  signal.  This  is  known  as  mul¬ 
tipath,  and  an  overview  of  its  effects  is  provided  in  Section  2.5  along  with  possible 
modelling  methods. 

These  distortions  can  be  modeled  as  a  linear  time-invariant  system  according 
to  (2.1),  which  maps  the  transmitted  signal  x  (t)  to  received  signal  y  (t)  using  impulse 
response  h  (t)  [14] . 


y{t) 


x  (A)  hit  —  A)  d\ 


(2.1) 


This  section  examines  the  baseband  communications  model  presented  in  Fig¬ 
ure  2.1.  The  transmitted  symbol  sequence  is  a  series  of  discrete  symbols  denoted 
by  x{n).  The  digital  to  analog  converter  (DAC)  block  converts  the  discrete  series 
x  [n)  to  a  continuous  time  (ci)  waveform  with  symbol  impulses  occurring  every  Tsym 
seconds.  The  transmitted  signal  is  then  acted  on  by  the  effective  channel  impulse 
response  h^jjit)  which  contains  the  pulse  shape  f^ct\t),  the  physical  channel  impulse 
response  h^JS (t) ,  and  a  matched  filter  in  the  receiver  /(ct)(— £),  described  in  (2.2). 
The  received  signal  y^{t)  is  then  sampled  by  an  oversampling  factor  of  O  times  per 
Tsym  [8]. 
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Equation  (2.2)  shows  the  form  of  h^j(t).  It  is  important  to  note  that  h^  {t) 
can  also  contain  multipath  effects  characterizing  the  different  paths  the  signal  can 
take  to  the  receiver.  Also  note  that  the  pulse  shapes  used  in  this  thesis  are  both  real 
and  symmetric,  so  that  the  pulse  shape  and  the  matched  filter  are  identical  to  one 
another. 


,  (2.2) 

where 


f(ct\t)  =  Pulse  shape  (and  matched  filter)  , 
h^phys(f)  =  Physical  channel  impulse  response. 


In  the  communications  model  presented  in  Figure  2.1,  the  N  transmitted  sym¬ 
bols  are  represented  by  a  series  of  discrete  terms,  x{n).  The  signal  after  matched 
Ritering  at  the  receiver  can  be  written  as 

N- 1 

y{ct)(t )  =  ^  x(n)h^f  (t  -  nTsym  -  tQ)  +  b{ct\t),  (2.3) 

n= 0 

where  is  the  effective  channel  impulse  response,  tQ  is  the  sampling  synchronization 
offset,  and  b^(t)  is  an  additive  noise  signal. 

The  additive  noise  signal  b^ct\t)  is  the  environmental  noise  through  the  trans¬ 
mission  channel  that  has  passed  through  the  matched  filter.  All  noise  used  in  this 
thesis  is  assumed  to  be  Additive  White  Gaussian  Noise  (AWGN),  denoted  as 
The  resulting  noise  signal  is  shown  in  (2.4). 


b^ct\t)  =  /  -  T)dT 


(2.4) 
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Due  to  the  oversampling  factor  O,  (2.3)  can  be  rewritten  as  (2.5)  with  t  replaced 
by  kTsym/ O  to  represent  the  higher  sampling  rate  of  Tsym/0  samples  per  second.  With 
this  substitution,  the  received  signal  after  matched  filtering  is  expressed  as 


The  discrete  analog  of  y ^  ^  fcT ^  is  given  in  (2.6)  and  represents  the  oversam¬ 
pled  discrete  received  sequence  before  equalization. 

! m  =  V{a)  (2.6) 

Equalization  attempts  to  minimize  or  remove  the  effects  of  ISI  and  other  channel 
effects  via  post  processing  within  the  communications  receiver.  This  can  be  done 
using  many  different  algorithms  and  techniques.  One  common  technique  used  for 
equalization  is  filtering.  Filtering  attempts  to  modify  the  received  signal  in  some 
way  that  makes  the  communications  symbols  distinct,  effectively  cancelling  out  or 
inverting  the  effect  of  the  channel  [14].  The  digital  transverse  filter  is  a  linear-time 
invariant  system  defined  by  the  discrete  impulse  response  w  which  acts  on  the  discrete 
received  symbols  y'{k)  to  produce  a  time  delayed  estimate  of  the  transmitted  symbols, 
defined  as  x  (k)  [11].  The  digital  transverse  filters’  operation  can  be  described  using 
the  discrete  convolution  sum, 

M—l 

x\k)  =  w'(k)  *y'(k )  =  ^  w'(m)y'(k  —  m)  .  (2.7) 

m= 0 

ffere,  x  (k)  is  the  oversampled  discrete  estimated  input  sequence  after  equaliza¬ 
tion,  and  is  then  downsampled  by  the  oversampling  factor  O  to  produce  the  estimated 
output  sequence  x(n).  Note  that  x(n)  is  approximately  equal  to  x{n  —  5),  and  5  is 


some  fixed  delay  that  maps  the  nth  estimated  symbol  to  the  nth  minus  5  transmitted 
symbol. 

The  mathematical  model  for  the  filtering  operation  given  in  (2.7)  is  only  valid  if 
the  channel  response  is  static,  and  in  most  real  physical  systems  the  channel  impulse 
response  varies  as  a  function  of  time  [14].  For  a  non  time  varying  channel,  a  preset 
equalization  method  is  sufficient  to  equalize  the  channel.  In  a  preset  equalization 
method,  the  equalization  taps  needed  to  invert  the  static  channel  are  calculated  and 
then  applied  to  each  received  signal,  with  no  need  to  adjust  the  filter  taps.  However, 
for  time  varying  channels,  some  method  of  adjusting  the  filter  taps  over  time  to 
account  for  changing  channel  conditions  is  needed.  An  adaptive  equalizer  is  one  that 
has  the  ability  to  update  the  filter  taps  either  continuously  or  at  set  periodic  times. 
The  convergence  rate  is  the  key  parameter  of  an  adaptive  filter  method  in  regards  to  its 
ability  to  equalize  a  time  varying  channel.  The  convergence  rate  is  a  measure  of  how 
quickly  the  adaptive  method  can  change  the  filter  taps,  and  if  the  channel  variation 
exceeds  the  convergence  rate  of  the  adaptive  filter,  symbol  estimation  performance  is 
degraded  [14]. 

A  subset  of  non-cooperative  receivers  using  blind  adaptive  equalization  tech¬ 
niques  are  examined  in  this  thesis.  These  techniques  are  designed  to  mitigate  the 
channel  effects  with  little  or  no  knowledge  of  the  channel  by  utilizing  the  statistics 
of  the  source  signal  [5].  These  types  of  blind  equalization  techniques  are  ideal  for 
use  in  the  design  of  a  non-cooperative  receiver  developed  for  eavesdropping,  and  the 
dither  algorithm  detailed  in  Section  2.6  varies  the  channel  impulse  response  beyond 
the  convergence  rate  of  these  blind  techniques  to  prevent  meaningful  third  party  de¬ 
modulation  of  the  transmitted  signal. 

2.3.1  Convolution  Matrix  Model.  The  discrete  linear  convolution  process 
described  in  Section  2.3  can  also  be  expressed  in  matrix  form.  As  convolution  is  a 
computationally  expensive  calculation,  recasting  the  convolution  process  as  a  series  of 
matrix  operations  allows  for  faster  and  more  efficient  computation.  It  is  important  to 
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Noise 

b^\n) 


note  that  the  baseband  model  presented  in  Figure  2.1  can  equivalently  be  represented 
as  a  multiple  channel  system,  with  each  of  O  subchannels  consisting  of  Tsym  spaced 
samples.  This  structure  is  shown  in  Figure  2.2.  For  modelling  a  communications 
system  as  a  series  of  matrix  computations,  the  effective  channel  impulse  response  can 
be  represented  as  a  convolutional  matrix,  and  can  aid  in  gaining  an  understanding 
of  the  equalization  process.  The  following  derivation  of  a  matrix-focused  model  is 
derived  from  [6]  and  [10]. 

Each  subchannel  impulse  response  h^l\n)  and  subchannel  equalization  filter 
w^\n)  sequence  can  be  defined  as  vectors  shown  in  (2.8)  and  (2.9),  where  L  is  the 
effective  channel  length,  M  is  the  number  of  equalizer  taps,  and  O  is  the  oversampling 
factor. 


h«  =  [hw( 0),  hw(  1),  ...  h{i\L  -  1)]T  ,  i  =  1, ...,  O  (2.8) 

WW  =  [ww(0),  ww(l),  ...  w®(M/0  -  1)]T  ,i  =  1, ...,  O  (2.9) 
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These  subchannel  equalization  vectors  can  then  be  combined  and  stacked  col¬ 
umn  wise  to  form  a  single  new  vector  w(n).  This  vector  is  formatted  differently  than 
the  oversampled  equalization  vector  w\  and  it  is  important  to  note  that  they  are  not 
equivalent. 


w(n) 


(2.10) 


The  matched  filter  output  at  the  receiver  of  each  subchannel  y^'in)  can  be 
expressed  as  a  vector  by  recasting  all  of  the  previous  discrete  sequences  as  vectors 
or  convolution  matrices.  Let  y W(n)  be  a  vector  of  N  samples  at  the  matched  filter 
output  as  shown  in  (2.11). 


yw(n)  =  [yM(n),  yM(n  -  1),  yM(n  -  2),  ...  y{l\n  —  N  +  1)]T  (2.11) 


The  transmitted  symbol  vector  x(n)  and  the  subchannel  AWGN  noise  vector, 
b«  (n),  are  defined  similarly  as  shown  in  (2.12)  and  (2.13). 


x(n)  =  [x(n),  x(n  —  1),  x(n  —  2),  ...  x(n  —  N  —  L  +  2)]T  ,  (2.12) 

b^(n)  =  [b^\n),  b^\n  —  1),  b^l\n  —  2),  ...  b^\n  —  N  +  1)]T  (2.13) 


The  subchannel  convolution  matrix  for  each  subchannel  is  a  Toeplitz  matrix 
with  the  subchannel  impulse  response  vector  h-*)  along  the  hrst  row  as  shown  in  (2.15). 
The  convolution  matrix  for  each  subchannel  acts  on  ( N  —  L  —  1)  transmitted  symbols 
in  x(n)  to  provide  N  samples  in  the  subchannel  output  vector  yW(n). 
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(2.14) 


y  ^(n)  =  H^x(n)  +  b^^n) 
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Stacking  all  subchannel  vectors  columnwise  allows  the  entire  multichannel  out¬ 
put  vector  y  (n)  to  be  expressed  in  terms  of  the  stacked  convolution  matrices  of  each 
subchannel: 


y(n)  = 


y(1)(n) 


y(°)(n) 


H 


(i) 

N 


H 


(O) 

N 


x  in) 


b«( 


n) 


b(°)( 


n) 


=  HArx(n)  +  b(n) 


(2.16) 


Similarly,  each  subchannel  equalization  filter  can  be  cast  as  a  convolution 
matrix  as  in  (2.15),  shown  in  (2.17).  Combining  these  convolutional  matrices  row¬ 
wise  as  shown  in  (2.18)  provides  a  method  to  calculate  an  estimate  of  the  transmitted 
vector  x(n),  shown  in  (2.19). 
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(2.17) 


vv  TV  vv  TV 


W  N  = 
x(n)  =  WNy{n 


(2.18) 

(2.19) 


2-4  Equalization  Techniques 

This  section  will  examine  the  both  the  cooperative  and  non-cooperative  equal¬ 
izer  algorithms  used  in  this  thesis,  and  provide  an  overview  of  the  important  charac¬ 
teristics  of  each.  An  overview  of  important  equalization  concepts  is  presented  first, 
covering  the  broad  types  of  equalizers  and  their  strengths  and  weaknesses.  The  coop¬ 
erative  receiver  structure  used  in  this  thesis  is  presented  next,  followed  by  an  overview 
of  non-cooperative  or  blind  equalizers  used. 

2-4-1  Equalization  Algorithm  Characteristics.  Blind  equalization  algorithms 
can  be  grouped  into  two  broad  categories:  algorithms  that  are  based  on  the  estimation 
of  a  transmitted  signal’s  second  order  statistics  (SOS),  and  algorithms  that  calculate 
equalizer  coefficients  based  on  higher  order  statistics  (HOS).  As  second  order  statis¬ 
tics  can  be  accurately  estimated  over  fewer  observations,  SOS  algorithms  typically 
converge  faster  than  algorithms  making  use  of  HOS  [5].  In  a  wireless  transmission 
environment  the  channel  conditions  often  vary  over  time,  making  algorithms  with  a 
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quick  convergence  rate  desirable.  If  the  channel  conditions  are  changing  more  rapidly 
than  an  equalizer  can  track,  then  the  equalizer  is  unable  to  invert  the  channel  to 
cancel  out  multipath  effects  making  data  extraction  difficult  [18]. 

Equalizers  can  also  be  classified  by  the  sampling  rate  used.  An  algorithm  typi¬ 
cally  either  uses  a  single  sample  from  each  communications  symbol  in  a  transmitted 
signal,  or  multiple  samples  from  each  symbol.  Algorithms  using  multiple  samples 
from  each  transmitted  symbol  are  known  as  fractionally  spaced  equalizers  (FSE), 
while  algorithms  using  a  single  sample  are  known  as  sample  or  baud  spaced  equal¬ 
izers  (BSE).  Fractionally  spaced  equalizer  algorithms  typically  perform  better  than 
equivalent  BSE  methods  [14],  converging  faster  and  more  accurately  [5]. 

Another  classification  parameter  worth  considering  is  an  algorithm’s  ability  to 
adapt  its  equalizer  coefficients.  As  stated  above,  the  wireless  transmission  environ¬ 
ment  often  includes  time  varying  channel  conditions.  Algorithms  that  update  the 
equalization  filter  after  every  received  symbol  are  known  as  continuously  adaptive. 
Algorithms  which  process  a  block  or  group  of  symbols  before  calculating  an  equal¬ 
ization  filter  are  known  as  periodically  adaptive.  An  algorithm  can  also  be  either 
a  direct  equalization  method  or  and  indirect  equalization  method.  A  direct  method 
calculates  the  filter  taps  for  a  linear  equalization  filter  which  is  then  applied,  where  an 
indirect  filter  estimates  the  channel  impulse  response  before  attempting  to  calculate 
an  equalization  filter. 

2.4.2  Minimum  Mean  Square  Error  (MMSE)  Equalization  .  The  coop¬ 
erative  equalizer  used  in  this  thesis  is  an  example  of  a  minimum  mean  squared  er¬ 
ror  (MMSE)  equalizer,  based  on  the  derivation  in  [12].  This  MMSE  equalizer  attempts 
to  minimize  the  cost  function  J(n)  shown  in  (2.20), 
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J(n)  =  E  [x*(n  —  5)x(n  —  <5)]  —  2Re  {e  \x(n  —  5)yi/(n)]  w*} 

+  wrE  [y(n)yH(n)]  w* 

=  E  [x*(n  -  8)x(n  -  5)]  -  2Re  {p^w*}  +  wT R?/w*  (2.20) 


where  w  is  the  equalization  filter,  Ry  is  the  autocorrelation  function  of  the  received 
vector,  and  pxy  is  the  cross  correlation  of  the  transmitted  and  received  signals.  Setting 
the  cost  function  gradient  given  in  (2.21)  to  zero  allows  for  the  optimal  mean  squared 
error  equalization  filter  to  be  calculated  exactly  according  to  (2.22). 


VwJ(n)  =  —2pxy  +  2wrRy 

W  opt  =  (Rj )  1  P  xy 


(2.21) 

(2.22) 


In  a  cooperative  receiver  context,  the  autocorrelation  and  cross  correlation  ma¬ 
trices,  R:y  and  pxy  respectively  are  easily  calculated,  and  the  implementation  of  this 
cooperative  method  in  this  thesis  in  the  presence  of  timing  dither  and  multipath  is 
presented  in  Chapter  3. 


2.4.3  Constant  Modulus  Algorithm  (CMA)  Equalization  .  The  constant 
modulus  algorithm  (CMA)  is  a  very  common  blind  equalization  technique  [6],  and 
both  baud  and  fractionally  spaced  versions  have  been  studied  and  documented  ex¬ 
tensively  in  literature.  CMA  uses  higher-order  statistics  to  directly  and  adaptively 
equalize  the  channel,  and  the  fractionally  spaced  version  provides  perfect  equaliza¬ 
tion  under  noiseless  conditions  for  most  channels  [6].  The  baud  spaced  version  of 
CMA,  indeed  any  BSE,  requires  an  infinite  equalizer  filter  length  to  achieve  perfect 
equalization,  making  it  somewhat  impractical  [4], 
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CMA  works  by  minimizing  an  error  term  e(n)  calculated  by  taking  the  difference 
between  the  expected  average  amplitude  and  the  sampled  magnitude  of  a  symbol,  as 
shown  in  (2.23). 


e{n)  = 


E[kW|4] 

E[|i(n)|2] 


x(n)  | 


(2.23) 


The  error  term  can  be  minimized  using  the  update  strategy  shown  in  (2.24), 
where  /i  is  the  convergence  step  size  of  the  algorithm,  w(n)  is  the  equalizer  coefficient 
vector,  and  y*(n)  is  a  regressed  vector  of  the  received  symbol  samples.  In  effect, 
the  CMA  algorithm  updates  its  nth  equalizer  filter  tap  by  operating  on  a  range  of 
received  symbols. 


w(n  +  1)  =  w(n)  +  /iy*(n)  (x(n)  —  yT(n)w(n))  , 

=  w(n)  +  /iy*(n)e(n)  .  (2.24) 

Being  a  minimization  acting  on  the  square  of  the  error  term,  e(n),  which  also 
contains  the  squared  symbol  estimate,  the  convergence  of  CMA  relies  on  the  estima¬ 
tion  of  fourth  order  statistics.  This  dependence  causes  CMA  to  have  a  low  convergence 
speed  when  compared  to  SOS  based  algorithms  [5].  If  the  signal  is  constant  modulus, 
meaning  that  each  symbol  has  the  same  amount  of  energy,  the  error  term  approaches 
zero  as  the  equalization  filter  converges.  The  CMA  algorithm  will  also  converge  for 
non  constant  modulus  signal  constellations,  but  will  converge  much  more  slowly  due 
to  the  cost  surface  having  a  non- zero  minima  [6]. 

CMA  assumes  a  wide  sense  stationary  (WSS)  or  cyclostationary  received  signal; 
WSS  for  a  baud  spaced  implementation  and  cyclostationary  for  a  FSE  implementa¬ 
tion.  Under  time  varying  channel  conditions  these  assumptions  will  be  violated,  but  if 
the  channel  impulse  response  varies  slowly  enough,  the  statistics  of  the  signal  may  still 
be  stable  enough  for  CMA  to  equalize  the  channel  [5].  In  FSE  CMA,  four  conditions 
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are  needed  to  guarantee  equalizer  filter  coefficients  that  perfectly  equalize  the  channel 
impulse  response.  First,  the  channel  must  be  noiseless.  The  source  signal  must  be  a 
series  of  independent  and  identical  symbols  selected  from  a  circularly  symmetric  con¬ 
stellation,  and  the  kurtosis  of  the  symbol  distribution  must  be  less  than  the  kurtosis 
of  a  Gaussian  distribution  [6].  Also,  the  channel  convolution  matrix  defined  by  Tong 
et  al.  in  [16]  must  be  full  column  rank. 

The  full  column  rank  condition  for  perfect  equalization  leads  to  a  minimum 
number  of  equalizer  filter  taps  required  in  implementing  the  CMA  algorithm.  The 
number  of  taps  needed  is  a  function  of  the  length  of  the  channel  response  and  the 
number  of  samples  taken  per  symbol,  or  the  sampling  rate,  and  is  shown  in  (2.25) 

M 

M>L  +  —  -1  (2.25) 

where  M  is  the  number  of  equalization  filter  taps,  L  is  the  effective  length  of  the 
channel  impulse  response,  and  O  is  the  oversampling  factor  [6]. 

2. 5  Multipath 

Multipath  interference  is  an  issue  that  commonly  arises  in  wireless  communica¬ 
tions.  A  large  body  of  research  exists  on  the  cancellation  and  equalization  of  multipath 
interference,  but  this  thesis  examines  the  effects  of  this  interference  on  an  ISI  time 
dithered  transmitter  and  receiver  pair.  This  ISI  waveform  necessitates  a  new  receiver 
structure  design  for  acceptable  demodulation.  This  section  begins  with  a  discussion 
of  the  causes  of  multipath  interference,  introduces  some  of  the  differing  effects  and 
systems  that  can  be  affected,  and  presents  some  possible  models  that  can  be  adapted 
to  the  ISI  multipath  problem  at  hand. 

As  the  name  implies,  multipath  interference  is  caused  by  having  multiple  time 
delayed  versions  of  a  transmitted  signal  arriving  at  the  receiver  as  a  result  of  the 
transmission  taking  several  different  paths.  This  can  affect  many  different  systems, 
but  radio  or  wireless  transmission  media  typically  have  the  worst  problem  with  mul- 
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tipath  due  to  the  difficulty  in  focusing  all  of  the  transmitted  energy  along  a  single 
path  [3].  Table  2.1  summarizes  the  mechanisms  involved  in  the  observed  multipath 
effects  of  several  different  systems. 


Table  2.1: 

Multipath  mechanisms  [3] 

System 

Multipath  mechanism 

HF  Radio 

Reflection  from  multiple  ionospheric  layers 

Mobile  and  personal  radio 

Reflection  and  scattering  from  buildings,  etc 

Microwave  point  to  point  links 

Atmospheric  refraction  and  reflection 

Satellite-mobile  systems 

Ground  and  building  reflection 

Radio  LAN/indoor  radio 

Reflections  from  walls  and  building  structure 

Diffuse  infrared 

Reflections  from  walls 

Multimode  optical  fibre 

Multimode  propogation 

Telephone/cable  network 

Reflections  from  terminations 

There  are  several  factors  influencing  the  amount  of  multipath  interference  present 
in  any  system,  and  the  steps  needed  to  attempt  to  equalize  it.  The  first  area  to  con¬ 
sider  is  the  characteristics  of  the  physical  transmission  channel.  The  geography  of 
the  transmission  area  may  hold  tall  buildings,  water  towers,  and  other  obstacles  that 
would  reflect  the  transmitted  signal  and  create  a  delayed  version  at  the  receiver  [2] .  If 
a  strong  line  of  sight  (LOS)  path  is  present  in  the  channel,  the  presence  of  such  a  path 
simplifies  the  equalization  problem,  but  in  some  environments  such  as  in  urban  areas 
the  LOS  path  is  often  blocked  or  obscured  by  surrounding  buildings.  If  the  trans¬ 
mitter  or  receiver  are  moving  relative  to  one  another,  then  any  multipath  reflections 
present  will  have  time-varying  delays  at  the  receiver.  These  time  delays  will  affect  the 
phase  of  the  transmitted  multipath  signal  with  respect  to  the  original  transmission, 
and  result  in  phase  distortion.  The  distance  between  the  transmitter  and  receiver 
plays  a  factor,  as  the  signal  suffers  free  space  path  loss  attenuation  according  to  the 
1/d2  law  [2],  Even  the  characteristics  of  the  objects  in  the  transmission  environments 
can  have  an  effect.  Some  objects,  such  a  a  metallic  water  tower,  will  reflect  more 
of  the  transmitted  energy  than  a  stone  building,  resulting  in  a  multipath  containing 
more  of  the  original  signal  power.  This  is  known  as  the  reflectivity  property  of  the 
object.  In  general,  a  transmission  environment  that  is  stationary,  meaning  that  the 


18 


transmitter /receiver  pair  and  objects  in  the  environment  are  not  moving,  is  easier  to 
equalize  at  the  receiver  than  time-varying  multipath  conditions. 

2.5.1  Multipath  Modelling.  In  order  to  develop  a  model  for  multipath  in¬ 
terference  for  use  in  simulations,  the  following  derivation  is  presented  based  upon  the 
development  in  [13].  First,  consider  a  transmitted  signal  s(t)  modulated  at  a  carrier 
frequency  fc  consisting  of  a  series  of  extremely  short  pulses.  s(t)  takes  the  form  shown 
in  (2.26) 


s(t)  =  $t[Sl(t)ej2*fct]  (2.26) 

where  si(t)  represents  the  baseband  transmitted  signal,  or  the  signal  before  modula¬ 
tion.  After  transmission,  multiple  propagation  paths  are  assumed,  and  each  path  has 
an  associated  attenuation  factor,  an{t)  and  delay  term,  rn{t).  Both  of  these  factors 
are  dependent  on  time  due  to  the  varying  nature  of  the  transmission  environment  [13]. 
The  received  bandpass  signal  at  the  receiver  x(t)  is  given  in  (2.27). 


x(t)  =  y;  an(t)s(t  -  rn(t ))  (2.27) 

n 

Substituting  s{t)  from  (2.26)  into  (2.27)  yields  the  new  complete  form  of  x{t) 


x[t)  =  9?< 


(2.28) 


representing  the  complete  model  of  the  received  bandpass  signal.  By  applying  lowpass 
filtering  to  x(t)  as  shown  in  (2.28),  the  lowpass  received  signal  r/(f)  can  be  written 
equivalently  as  (2.29). 


nit)  =  Y  an{t)e-^^Sl{t  -  rn(t ))  (2.29) 

n 


19 


Since  ri(t)  is  the  response  of  a  lowpass  channel  to  the  equivalent  lowpass  signal, 
si(t),  the  lowpass  channel  can  be  described  by  the  time  variant  impulse  response  given 
in  (2.30)  [13].  This  also  represents  the  multipath  channel,  (t),  as  given  in  (2.2). 


c(t;  t)  =  ^2an(t)e  -  Tn(t)) 


(2.30) 


n 


Assuming  the  lowpass  channel  impulse  response  c(r;  t)  is  wide  sense  station¬ 
ary  (WSS),  the  the  autocorrelation  function  of  c(r;  t)  is  defined  as  (j)c. 


0c(T1}T2;Ai)  =  i£'[c*(r1;t)c(r2|(  +  At)] 


(2.31) 


As  each  propagation  path  passes  through  a  different  transmission  environment, 
each  is  subject  to  different  attenuations  and  phase  shifts.  In  most  wireless  transmis¬ 
sion  media,  these  attenuations  and  phase  shifts  are  uncorrelated  with  one  another  [13]. 
Assuming  these  factors  are  uncorrelated  for  two  different  delays,  <pc  can  be  written  as 
shown  in  (2.32). 


[c*(ti;  t)c(r2;  t  +  At)]  =  0c(ri;  At)S(ri  -  t2) 


(2.32) 


Letting  At  =  0,  the  autocorrelation  function  </>c(r;  0),  or  0c(r),  is  the  average 
power  output  of  the  channel  as  a  function  of  time.  As  such,  </>c(r)  is  known  as  the 
multipath  intensity  profile  of  the  channel  [13].  In  many  typical  systems,  the  measured 
4>c(t)  appears  as  shown  in  Figure  2.3.  Physically,  this  represents  the  expected  power 
of  a  multipath  coefficent  of  a  channel  at  a  particular  delay  r. 

2.6  Dither  Algorithm 

This  section  provides  an  overview  of  the  method  developed  in  [8]  to  induce 
time  varying  ISI  into  a  transmitted  waveform  to  make  it  difficult  to  equalize.  As  the 
equalization  techniques  discussed  in  Section  2.4  all  rely  on  a  slowly  changing  effective 
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Figure  2.3:  Typical  Multipath  Intensity  Profile 


channel  impulse  response,  any  ISI  injected  will  have  to  be  severe  enough  to  prevent 
third  party  demodulation.  The  time  varying  aspect  of  the  induced  ISI  must  be  large 
enough  to  frustrate  blind  equalization  techniques  as  well,  but  not  so  severe  as  to 
prevent  effective  cooperative  equalization. 

Most  traditional  receivers  transmit  one  symbol  at  a  time,  with  each  new  symbol 
beginning  at  some  set  time  spacing,  Tsym,  as  stated  in  Section  2.3.  To  create  the 
ISI  in  the  waveform,  the  beginning  of  each  symbol  is  delayed  within  the  range  of 
TSym-  The  dither  waveform  can  then  be  pictured  as  a  series  of  Tsym  time  blocks  which 
each  contain  the  leading  edge  of  exactly  one  communications  symbol.  This  limitation 
preserves  the  traditional  transmission  symbol  rate,  and  simplifies  the  design  of  a 
cooperative  receiver  [8]. 

This  delay  is  applied  to  the  transmitted  waveform  immediately  after  the  DAC, 
as  shown  in  Figure  2.4,  with  the  delay  of  a  particular  symbol  n  from  the  start  of  the 
nth  Tsym  and  is  defined  as  d(n).  After  the  application  of  this  dither,  the  transmitted 
signal  is  expressed  in  (2.33) 

AT— 1 

=  ^2  x(n)S  (t  -  nTsym  -  d{n))  ,  (2.33) 

n= 0 

where  x(ri)  is  the  transmitted  data  symbols,  and  d(n)  is  a  sample  of  a  random  vari¬ 
able  D.  D  is  a  uniform  random  variable  with  probability  density  function  of  p^\t). 
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xyn) 
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Figure  2.4:  A  dither  block  in  the  transmitter  creates  ISI  [8] 


The  range,  Tmax,  of  the  uniform  random  variable  D  can  be  selected,  so  long  as  it 
falls  within  [0,  Tsym\.  The  probability  density  function  can  be  written  in  terms  of  the 
selected  range  Tmax  as  shown  in  (2.34). 


I  I 


0, 


<  0  or  t  >Tr 


(2.34) 


To  create  a  sequence  of  pseudorandom  delays  d(n),  each  individual  delay  must 
be  selected  from  a  discrete  set  of  possible  delays.  Based  on  a  uniform  distribution, 
the  symbol  period  is  divided  into  P  +  1  equally  spaced  possible  delay  values.  The 
discrete  probability  density  function  Po(t)  is  given  in  (2.35) 


PD{t) 


Pn 


+  1 


Pmax 

T.s 

p=0 


(2.35) 


where  pmax  is  the  largest  integer  delay  value  contained  within  Tmax.  The  relationship 
between  pmax  and  Tmax  is  given  in  (2.36),  where  [J  is  the  floor  function  which  rounds 
down  to  the  largest  integer  value. 


Pmax 


PT 

1  J-rr 


■  sym 


(2.36) 


As  a  finite  sequence,  the  delay  code  d{n)  will  have  some  associated  period,  N^. 
In  all  cases,  this  thesis  considers  long  delay  codes  with  periods  larger  than  the  number 
of  transmitted  data  symbols.  It  is  also  possible  to  develop  cooperative  receivers  that 
can  function  provided  a  short  period  code  is  used  relative  to  the  number  of  transmitted 
symbols  available  for  equalizer  training  [8],  but  the  performance  of  such  a  receiver  in 
the  presence  of  multipath  is  not  examined  in  this  thesis.  Excluding  the  effects  of 
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multipath  interference,  the  severity  of  the  ISI  created  through  the  dither  process  is  a 
function  of  the  pulse  shape  characteristics  and  the  range  of  Tsym  over  which  the  delay 
probability  density  function  is  allowed  to  use,  Tmax.  In  this  thesis  various  fixed  pulse 
shapes  will  be  used  with  varying  symbol  dither  to  characterize  the  non-cooperative 
suppression  capabilities  of  the  developed  dither  algorithm  in  the  presence  of  multipath 
interference. 
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III.  Research  Methodology 


3.1  Overview 

This  chapter  presents  a  numerical  simulation  method  for  testing  the  performance 
of  the  ISI  algorithm  reviewed  in  Section  2.6,  in  the  presence  of  multipath  in¬ 
terference.  Section  3.2  introduces  the  overall  model  used  to  build  up  the  simulation 
designed  and  implemented  in  MATLAB,  with  Sections  3.3,  3.4,  and  3.5  discussing 
the  transmitter  model,  the  channel  model,  and  the  receiver  model  respectively.  The 
binary  DPSK  constellation  used  in  this  thesis  is  detailed  in  Section  3.3,  along  with 
the  pulse  shapes  used  for  the  communications  symbols.  A  derivation  of  the  multipath 
interference  models  is  provided  in  Section  3.4,  and  Section  3.5  includes  derivations 
of  how  the  equalization  models  reviewed  in  Chapter  11  are  implemented  in  simula¬ 
tion.  Section  3.6  provides  a  progressive  validation  of  the  different  components  of  the 
numerical  simulation  model. 


3.2  Model  Overview 


The  end-to-end  complete  communications  system  modeled  in  Matlab  in  this 
thesis  is  shown  in  Figure  3.1.  The  transmitter  block  takes  the  input  bit  stream 
and  creates  a  baseband  waveform  of  binary  DPSK  symbols.  The  channel  block  injects 
AWGN  and  multipath  interference  effects  into  the  waveform,  and  the  receiver  block 
estimates  the  transmitted  symbols  to  produce  the  output  bit  sequence  bout.  As  the 
intent  of  this  thesis  is  to  examine  the  performance  of  the  dither  algorithm  developed 
in  [8],  the  original  transmitter  structure  remains  unchanged.  The  channel  model 
and  correspondingly  the  receiver  models  have  been  adapted  to  include  multipath 
interference  explicitly  as  detailed  in  the  corresponding  sections  below. 


{bin} 


Figure  3.1:  Simulation  Block  Diagram  [8] 
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3.3  Transmitter  Model 


The  transmitter  block  diagram  shown  in  Figure  3.2  details  the  steps  involved  in 
creating  the  sampled  transmitted  waveform,  u(p).  The  simulation  uses  binary  DPSK 
modulation  to  encode  the  input  bit  sequence  into  a  series  of  symbols. 


x(n)  xup(p) 

Figure  3.2:  Simulated  transmitter  block  diagram  [8] 


The  DPSK  modulation  is  circularly  symmetric  about  the  origin  in  the  signal  or  con¬ 
stellation  space,  with  equal  magnitude  symbols  and  uniform  phase  separation  between 
symbols.  In  DPSK,  a  constellation  of  M  symbols  is  separated  into  a  series  of  words  of 
length  /,  where  l—log^M)  bits.  These  l  bit  words  are  mapped  to  a  phase  offset  using 
Gray  coding.  Figure  3.3  shows  the  binary  DPSK  constellation,  and  Table  3.1  provides 
the  corresponding  Gray  code  symbol  mapping  [14].  Binary  DPSK  also  satisfies  the 
conditions  for  perfect  equalization  of  the  CMA  algorithm  as  discussed  in  Chapter  II. 
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Figure  3.3:  Binary  DPSK  Constellation 


Real 


The  upsampling  block  increases  the  sampling  frequency  by  oversampling  by  a 
factor  of  P,  (P=64  throughout  this  thesis)  placing  an  integer  number  of  zero  valued 
samples  between  each  transmitted  symbol.  The  signal  is  now  a  sequence  of  symbols 
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Table  3.1:  Binary  DPSK  Gray  code  mapping 


Binary 

Word 

Phase  Shift,  A  9 

0 

0 

1 

7 r 

spaced  P  samples  apart.  An  upsampling  factor  of  64  is  used  throughout  this  thesis. 
Upsampled  signals  are  indexed  using  the  parameter  p.  Two  variables  are  used  for 
upsampling  and  downsampling  in  this  thesis:  P  and  O.  These  two  different  variables 
are  used  to  model  a  continuous  time  signal  using  a  discrete  time  model.  Upsampling 
by  P  ensures  that  the  continuous  time  signal  is  sampled  at  a  high  enough  rate  to 
approximate  the  signal  digitally.  The  downsampling  factor  0  is  used  in  the  receiver 
to  downsample  the  signal  by  32  to  make  use  of  Tsym/ 2  fractionally  spaced  equalizers. 
This  downsampling  factor  could  be  any  value  equal  or  less  than  64,  but  fractionally 
spaced  equalizers  provide  a  good  tradeoff  in  saved  complexity  and  equalizer  perfor¬ 
mance.  The  delay  block  implements  the  symbol  timing  dither  algorithm  overviewed 
in  Section  2.6,  moving  each  impulse  an  integer  number  of  samples  as  specified  by  the 
delay  code.  To  implement  the  symbol  dither,  the  dither  algorithm  generates  a  pseudo¬ 
random  delay  code  d  consisting  of  integer  delays  to  be  applied  to  each  corresponding 
symbol.  The  maximum  allowable  delay  must  be  less  than  P-1  samples  to  prevent 
transmitting  symbols  out  of  order.  The  distribution  of  delays  is  uniform  across  the 
maximum  allowable  delay  period  throughout  this  thesis,  and  this  maximum  delay 
term  can  be  adjusted  to  limit  the  amount  of  ISI  induced  in  the  transmitted  wave¬ 
form.  The  signal  is  then  convolved  with  the  pulse  shape  signal  f(p)  to  produce  the 
transmitted  waveform  u(jp). 

3.3.1  Pulse  Shape.  This  thesis  uses  a  Square  Root  Raised  Cosine  (SRRC) 
pulse  shape,  with  several  different  roll  off  factors  (3 .  The  /3  value  measures  the  amount 
of  energy  that  is  in  the  additional  lobes  of  the  SRRC  shape,  with  lower  beta  values 
having  more  energy  outside  the  main  lobe.  Figure  3.4,  shows  the  three  SRRC  shapes 
used  in  this  thesis.  As  the  SRRC  pulse  shape  is  unbounded,  it  must  be  truncated  in 
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any  physical  system.  Throughout  this  thesis,  the  SRRC  pulse  width  pw  is  10  symbol 
periods  wide  centered  at  t=0. 


Figure  3.4:  SRRC  pulseshapes  with  /3— 0,  0.25  &  0.5 
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3.4  Channel  Model 


The  simulated  physical  channel  used  in  this  thesis  is  a  combination  of  the  mul¬ 
tipath  conditions  present  in  the  channel  ( hphys )  and  an  AWGN  noise  signal  v(p)  as 
shown  in  Figure  3.5. 


u(p) 


Noise 

u(p) 


Physical 

Channel 


*■  r(p) 


Figure  3.5:  Simulated  channel  block  diagram  [8] 


After  applying  the  multipath  effects  of  the  channel,  the  AWGN  signal  is  added. 
The  noise  samples  are  complex,  with  the  real  and  imaginary  components  of  each 
sample  having  independent  Gaussian  distributions.  The  noise  sample  variances  are 
then  scaled  to  achieve  a  desired  SNR  at  the  input  of  the  receiver,  and  are  expressed 
in  terms  of  Eb/Na.  This  simulated  SNR  is  calculated  using  the  mean  squares  of  the 
transmitted  signal  z(p)  and  the  noise  signal  v(p)  as  shown  in  (3.1) 


SNR 


EKg)  I2 

E  Kp)  I2 


(3.1) 


It  is  also  important  to  note  that  in  the  presence  of  ISI  due  to  symbol  timing  dither, 
pulse  shapes,  and  multipath  conditions  can  be  treated  as  a  noise  source.  These  are 
not  included  in  SNR  calculations  in  this  thesis.  The  multipath  interference  models 
used  and  their  implications  are  discussed  in  Section  3.4.1. 


3-4-1  Multipath  Model.  The  purpose  of  this  thesis  is  to  examine  the  perfor¬ 
mance  of  a  symbol  timing  dither  waveform  in  the  presence  of  multipath  interference. 
In  order  to  do  that,  two  separate  multipath  profile  models  are  presented  here.  The 
first  is  composed  of  a  line-of-sight  path  modeled  with  zero  delay  and  one  multipath 
reflection  delayed  3 Tsym,  the  second  consists  of  a  line-of-sight  path  with  zero  delay 
and  3  multipath  reflections,  at  2 Tsym,  4 Tsym  and  6Tsym  delays,  respectively.  In  order 
to  provide  diversity  in  the  multipath  modelling,  50  distinct  channels  are  generated 
for  both  the  2  and  4  ray  models.  Each  of  these  channels  is  randomly  generated  by 
defining  a  deterministic  delay  profile  based  on  (3.2) 

m  =  e~at  (3.2) 

where  a  represents  a  scalar  modeling  the  rate  of  exponential  decay  [7].  For  the  2  ray 
model,  the  line-of-sight  channel’s  amplitude  is  defined  as  yielding  unity  power,  and 
the  multipath  reflection  power  is  V0.35.  For  the  4  ray  model,  the  three  multipath 
reflections  have  power  values  of  ^0.22,  \/0.05,  \/0.02,  respectively.  These  choices 
are  consistent  with  the  typical  multipath  intensity  profile,  shown  in  Figure  2.3.  Each 
multipath  reflection  power  value  is  then  multiplied  by  a  zero  mean  unit  variance  Gaus¬ 
sian  random  variable.  This  generates  50  unique  channels  for  each  multipath  model, 
where  each  multipath  reflection  is  normally  distributed  with  the  average  power  value 
selected  above.  The  delays  between  multipath  reflections  is  held  constant  for  all  sim¬ 
ulations  of  a  given  model  throughout  this  thesis.  Each  channel  is  then  normalized  to 
unity  power  to  prevent  adding  any  artificial  power  to  the  system  through  the  multi- 
path  interference.  To  produce  results  reflecting  the  average  performance  of  a  symbol 
timing  dither  waveform  in  multipath,  the  simulation  transmits  iV=1000  symbols  for 
training  the  CMA  equalizer,  and  1V=1000  symbols  for  estimation  using  each  of  these 
random  channels  in  turn  as  hphys-  The  bit  error  rate  (BER),  Pft  is  calculated  over  each 
channel  and  averaged  together  to  provide  the  final  results,  presented  in  Chapter  IV. 
Definitions  for  P&  used  in  this  thesis  are  provided  in  Section  3.5. 
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Figure  3.6:  Simulated  receiver  block  diagram  [8] 


3.5  Receiver  Model 

The  simulated  receiver  block  diagram  is  shown  in  Figure  3.6.  The  receiver  takes 
the  input  received  signal  r(p)  and  convolves  it  with  the  same  pulse  shape  used  in  the 
transmitter  block.  This  applies  the  effects  of  matched  filtering  at  the  receiver.  The 
resulting  signal  is  then  downsampled  by  a  factor  of  2/P.  This  results  in  having  2 
samples  from  every  transmitted  symbol  period,  enabling  fractionally  spaced  equal¬ 
ization  using  two  subchannels  in  both  the  cooperative  MMSE  and  non-cooperative 
CMA  algorithms  presented  in  this  thesis.  The  baseline  symbol  spaced  receiver  only 
uses  data  from  the  first  subchannel,  which  is  analogous  to  symbol  spaced  sampling. 
This  simulation  synchronizes  the  first  subchannel  to  the  point  of  the  maximum  effec¬ 
tive  channel  impulse  response  in  the  absence  of  symbol  timing  dither.  This  provides 
the  highest  SNR  values  for  the  basic  symbol  spaced  receiver.  In  the  presence  of  sym¬ 
bol  timing  dither,  the  simulation  synchronizes  to  the  point  of  the  average  maximum 
effective  channel  response  within  the  symbol  period. 

The  binary  DPSK  demodulator  uses  the  phase  difference  between  the  current 
symbol  and  the  previous  estimated  symbol  to  estimate  transmitted  information  bits. 
Because  this  process  is  dependent  only  on  the  phase  difference  between  adjacent 
estimated  symbols  and  not  the  absolute  phase  values,  this  technique  is  immune  to 
rotation  of  the  constellation  in  the  signal  space  [14].  For  DPSK,  the  probability  of  bit 
error  is 

A  =  (3.3) 
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where  E b  is  the  average  energy  per  bit  in  the  information  signal  and  Na  is  the  noise 
power  spectral  density.  All  analytical  Pb  curves  presented  in  this  thesis  are  generated 
using  (3.3). 

3.5.1  Cooperative  Receiver.  In  a  known  channel,  MSE  minimization  is  a 
common  approach  for  calculating  an  equalization  filter  to  deconvolve  channel  effects. 
Given  a  constant  or  static  channel  impulse  response,  the  MSE  minimization  process 
produces  one  single  equalization  filter.  As  the  multipath  interference  present  in  a 
channel  is  held  constant  for  the  purposes  of  equalization  in  this  simulation,  the  mul¬ 
tipath  aspect  of  the  channel  is  stable.  However,  in  the  presence  of  symbol  timing 
dither,  each  element  of  the  delay  code  will  constitute  a  unique  channel  to  equalize, 
necessitating  the  computation  of  Nj  equalizers. 

The  basic  structure  of  the  MSE  algorithm  is  presented  in  Section  2.4.2,  with 
the  optimal  MSE  equalization  filter  given  by  (2.22).  In  adapting  the  MSE  algorithm 
for  use  in  a  pseudorandom  delay  waveform  with  multipath,  a  fractionally  spaced 
implementation  is  used.  To  signify  the  ISI  caused  by  the  symbol  timing  dither,  let 
H  N  (n)  represent  the  convolution  matrix  that  will  map  the  transmitted  symbols,  x(n) 
to  the  output  vector  y(n),  similar  to  the  formation  shown  in  (2.16).  As  each  symbol 
can  potentially  have  a  unique  delay  due  to  the  symbol  timing  dither,  the  convolution 
matrix  must  also  be  indexed  by  n. 

The  autocorrelation  of  the  output  signal  is  given  in  (3.4),  and  is  also  dependent 
upon  the  symbol  timing  dither  applied  to  each  symbol  and  is  indexed  by  n. 


R  y(n) 


=  E  [y(n)yH(ra)]  , 


=  E 


(Hjy(n)x(n)  +  b(n))  (Hjv(n)x(n)  +  b(n)) 


H 


=  E  [Hyv(n)x(n)x^(n)H^(n)]  +  E  [b^b^ra)]  , 


=  Hjv(n)R;cH^(n)  +  Rb  . 


(3.4) 
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Again,  the  cross-correlation  vector  p xy(n)  is  dependent  on  the  symbol  index  n: 

pxy(n)  =  E[y*(n)x{n- 5)}  , 

=  E  [(Hjv(^)x(n)  +  b(n))*  x(n  —  5)]  , 

=  H *N(n)E  [x*(n)x(n  —  5)}  .  (3.5) 

Finally,  the  MSE  optimal  filter  must  now  indexed  as  well: 

wopt(n)  =  (Rj'(n))  1p*y(n)  .  (3.6) 

Equation  (3.6)  shows  the  form  of  the  optimal  MSE  filter  that  must  be  calculated 
for  each  symbol.  In  order  to  calculate  the  optimal  filter,  the  terms  defining  Ry(n) 
and  p xy(n),  namely  Hjv(u),  R^R*,,  and  E  \x*(n)x{n  —  5)]  must  be  defined.  Based  on 
the  derivations  found  in  Section  3.2.2  of  [8],  these  parameters  are  presented  below. 


H$(n)  = 


N+L 


h(n\  0) 


h^_L+1(L-l) 


o  e.fo) 


hn-L(L-  1) 


ln-N+l 


(0) 


0 

0 

i  ( i ) 

'ln-N-L+ 2 


(L-  1) 


N 


(3.7) 


Equation  (3.7)  represents  the  subchannel  convolution  matrix  for  a  particular 
symbol  n.  Each  row  contained  in  H^- (n)  is  constructed  from  heff  as  given  by  2.2, 
where  hphys  includes  the  effect  of  the  selected  multipath  interference  channel.  These 
rows  also  account  for  the  ISI  generated  by  nearby  symbols  on  the  symbol  of  interest  n 
due  to  symbol  timing  dither.  Each  pulse  shape  is  10  symbol  periods  wide,  so  a  window 
of  n- 5  to  n+5  symbols  could  potentially  influence  the  nth  symbol.  The  convolution 
matrix  Htv(u)  is  constructed  by  stacking  the  subchannel  matrices  for  each  symbol 
columnwise. 
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Assuming  the  source  sequence  x(n)  is  zero  mean  and  ideally  and  identically 
distributed,  Rx  is  an  identity  matrix  scaled  by  E  [x*{n)x(n)\  and  E  [x*(n)x(n  —  <5)] 
has  precisely  one  nonzero  element  valued  as  E  [x*(n)x(n)}  at  the  (5  +  1)  element. 

The  elements  discrete  sampled  elements  of  can  be  defined  as  in  (3.8) 

E  [b'ik^b'*^)]  =  No9W  [(n2  -  m)Ts,m  +  (i2  -  h )Tsym/0]  (3.8) 

where  g^(t)  represents  the  convolution  of  the  pulse  shape  and  multipath  effects  and 
the  corresponding  matched  filter  scaled  by  the  noise  power  density,  Na. 


—  f(ct) 


( Ct)*( 


f(d)  (t)  {t  +  T)dT 


(3.9) 


These  parameters  can  then  be  used  with  (3.6)  to  calculate  each  of  the  neces¬ 
sary  equalization  filters  for  a  particular  delay  code.  As  long  as  the  multipath  profile 
and  the  delay  code  remain  the  same,  these  Nj  equalization  filters  can  be  reused  assum¬ 
ing  synchronization  is  maintained.  The  calculation  of  each  optimal  filter  requires  that 
the  cooperative  system  have  full  knowledge  of  the  pseudorandom  delay  sequence  used 
and  the  multipath  conditions  present  in  the  channel.  As  the  multipath  conditions  in 
the  channel  can  be  measured  before  actual  transmission  using  a  training  signal,  this 
is  feasible.  Also,  the  equalizer  must  have  some  set  delay  built  into  the  process,  i.e. 
the  estimated  symbol  x(n  —  5)  is  produced  by  the  corresponding  filter  w opt(n). 

3.5.2  Non-Cooperative  Receivers.  There  are  two  non-cooperative  receiver 
structures  presented  in  this  thesis:  a  baseline  symbol  spaced  receiver  and  a  frac¬ 
tionally  spaced  CMA  equalizer.  Both  are  assumed  to  have  knowledge  of  the  symbol 
constellation,  symbol  period  and  the  transmitted  pulse  shape,  but  have  no  knowledge 
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of  the  pseudorandom  symbol  timing  delay  values  provided  in  d(n).  The  performance 
of  both  equalizers  in  the  presence  of  multipath  and  symbol  timing  dither  is  used  to 
compare  performance  with  the  MMSE  cooperative  equalizer. 

The  basic  receiver  samples  the  matched  filter  output  at  the  receiver  at  a  rate 
of  1/Tsym.  This  is  analogous  to  symbol  or  baud  spaced  sampling.  The  bit  error 
rate  results  of  this  receiver  form  a  baseline  for  comparison  of  the  other  equalizer 
structures  in  the  presence  of  multipath  and  symbol  timing  dither;  as  no  equalization 
is  performed,  both  the  cooperative  and  non-cooperative  equalization  methods  should 
best  this  baseline  performance. 

The  FSE  CMA  equalizer  utilizes  the  Constant  Modulus  Algorithm  with  a  frac¬ 
tionally  spaced  implementation,  i.e. ,  the  matched  filter  is  sampled  at  2 /Tsym.  The 
CMA  equalizer  sampling  phase  is  also  synchronized  as  in  the  basic  receiver.  CMA 
equalization  is  detailed  in  Section  2.4.3.  The  CMA  algorithm  requires  that  the  over¬ 
sampled  matched  filter  outputs  be  cyclostationary  for  perfect  equalization,  which  is 
invalidated  through  the  ISI  induced  by  the  presence  of  symbol  timing  dither.  How¬ 
ever,  the  adaptive  nature  of  the  CMA  algorithm  could  still  lead  to  improved  bit  error 
rate  performance  through  equalizing  the  average  ISI  due  to  multipath  and  symbol 
timing  dither  conditions. 


34 


3. 6  Simulation  Validation 


This  section  presents  the  validation  process  results.  Each  figure  presented  in 
this  section  is  generated  using  the  deterministic  base  multipath  models,  and  include 
no  randomized  multipath  channels  as  described  in  Section  3.4.1.  Each  data  point 
is  generated  using  50,000  symbols  for  estimation.  All  confidence  intervals  shown  in 
the  validation  curves  are  95%  confidence  intervals  based  on  the  number  of  bit  errors, 
approximating  Monte  Carlo  simulation  methods. 

3.6.1  Assumptions.  This  section  points  out  the  assumptions  and  simplifica¬ 
tions  present  in  this  simulation.  Doppler  effects  are  not  modeled  in  this  simulation, 
as  the  frequency  of  the  carrier  is  assumed  to  be  equal  for  each  multipath  reflection. 

The  receiver  is  capable  of  finding  the  point  at  which  the  combined  channel 
impulse  response  is  at  a  maximum  and  synchronizes  to  sample  the  waveform  at  this 
point.  This  synchronization  is  critical  to  the  performance  of  the  basic  symbol  spaced 
receiver  structure.  After  the  pseudorandom  symbol  timing  dither  is  included  in  the 
system,  the  point  at  which  this  maximum  occurs  changes  from  symbol  to  symbol.  The 
receiver  will  then  synchronize  to  the  average  time  in  each  symbol  period  at  which  the 
impulse  response  is  at  a  maximum.  This  synchronization  effect  is  applied  within  the 
simulation,  and  is  assumed  to  be  ideal. 

The  cooperative  MMSE  equalizer  structure  requires  knowledge  of  the  AWGN 
noise  power.  In  simulation  this  is  determined  directly  using  the  variance  of  the  noise 
signal.  In  practice,  the  noise  signal  would  not  be  available  separate  from  the  trans¬ 
mitted  signal  at  the  receiver,  but  this  could  be  achieved  in  practice  by  determining 
the  noise  power  at  the  receiver  before  transmission  begins. 
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(dB) 

Figure  3.7:  Validation  with  no  timing  dither  and  SRRC  /3=0.25 

Figure  3.7  represents  a  baseline  simulation  to  test  model  validity.  No  multipath 
or  symbol  timing  dither  effects  are  included,  and  the  analytical  P&  is  derived  from  the 
ideal  binary  DPSK  performance  curves  as  defined  in  [14].  All  equalization  methods 
are  statistically  identical  to  the  analytical  curve  under  these  conditions,  as  expected 
for  a  valid  simulation  model. 

Figure  3.8  introduces  symbol  timing  dither  to  the  system,  with  a  max  timing 
delay  of  0.6 Tsym.  The  symbol  timing  dither  will  create  ISI  in  the  system,  degrading 
the  performance  of  all  equalizers.  As  the  cooperative  MMSE  equalizer  uses  knowledge 
of  the  pseudorandom  delay  code  generating  the  symbol  timing  delays,  it  performs 
better  in  terms  of  bit  error  rate  than  both  the  basic  symbol  spaced  receiver  and  the 
non-cooperative  CMA  algorithm  equalizer. 
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Eb/No  (dB) 

Figure  3.8:  Validation  with  0.6 Tsym  max  delay  and  SRRC  /3=0.25 


E/K  (dB) 


Figure  3.9:  Validation  with  two  ray  multipath  and  SRRC  /3=0.25 

Figure  3.9  and  Figure  3.10  show  the  simulation  performance  for  the  2  and  4 
ray  multipath  models  in  the  absence  of  symbol  dither,  respectively.  The  bit  error 
rate  performance  of  both  the  CMA  and  cooperative  MMSE  equalizers  is  in  the  region 
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between  the  symbol  spaced  receiver  performance  and  the  analytical  binary  DPSK  per¬ 
formance  curve,  indicating  that  both  equalizers  are  partially  equalizing  the  multipath 
interference  present.  The  CMA  algorithm  uses  an  adaptive  update  strategy  presented 
in  Section  2.4.3,  and  as  such  has  a  certain  amount  of  variation  in  the  equalizer  filter 
taps  derived  from  the  step  size,  /i.  This  could  be  a  contributing  factor  to  the  degraded 
Pb  performance  relative  to  the  cooperative  MMSE  equalizer  structure. 


H/N,  (dB) 

Figure  3.10:  Validation  with  four  ray  multipath  and  SRRC  /3= 0.25 

Figure  3.11  and  Figure  3.12  show  the  validation  results  for  the  2  and  4  ray 
models  with  symbol  timing  dither  effects  included.  In  both  of  these  figures  the  co¬ 
operative  MMSE  equalizer  structure  performs  better  than  its  non-cooperative  CMA 
counterpart,  providing  a  performance  improvement  of  approximately  4  dB  in  terms  of 
bit  error  rate  at  Pb  =  4  x  ICR2.  However,  the  MSE  equalizer  performance  also  suffers 
statistically  significant  degradation  from  the  analytical  Pb  curve.  These  validation 
results  are  consistent  with  prior  results  presented  in  [8],  and  indicate  that  the  various 
parts  of  the  simulation  model  are  functioning  correctly.  The  effects  on  bit  error  rate 
for  a  full  range  of  random  multipath  channels,  max  symbol  timing  delays,  and  pulse 
shapes  is  presented  in  Chapter  IV. 
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Figure  3.11:  Validation  with  two  ray  multipath,  0.6 Tsym  max  delay  and  SRRC 

(3=0.25 


H/N,  (dB) 


Figure  3.12:  Validation  with  four  ray  multipath,  0.6 Tsym  max  delay  and  SRRC 

/3=0.25 
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3.6.2  Simulation  Parameters.  This  section  provides  two  tables  summarizing 
the  simulation  parameters  used  in  this  thesis  for  the  reader’s  convienience.  Table  3.2 
provides  a  list  of  all  of  the  global  simulation  parameters,  while  Table  3.3  details  the 
parameters  used  by  both  the  cooperative  and  non-cooperative  receiver  structures  used 
in  simulation. 


Table  3.2:  Simulation  Parameters 


p 

Number  of  samples  per  symbol  period 

M 

Number  of  symbols  in  the  DPSK  constellation 

Pulse  Shape 

Pulse  shape  in  the  transmitter  and  the  corresponding 
matched  filter 

d 

Delay  code  sequence  used  to  determine  the  delay  of  each 
transmitted  symbol 

Vmax 

Maximum  allowable  delay  in  the  sequence  d 

Vi 

LLphys 

Physical  channel  model 

SNR 

Signal-to-noise  ratio  as  measured  the  receiver  input 

0 

Number  of  samples  per  symbol  period  after  receiver 
downsampling 
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Table  3.3:  Equalization  Parameters 


Baseline 

N/A 

No  parameters  necessary 

FS  CMA 

h 

Gradient  descent  step  size 

M 

Number  of  taps  in  the  equalizer  (equalizer  length) 

N 

Number  of  symbols  used  to  train  the  equalizer 

Modified  MMSE 

S 

Impulse  response  delay  to  which  the  equalizer  converges 
(in  symbol  periods) 

L 

Estimate  of  the  length  of  the  effective  channel  (in  symbol 
periods) 
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IV.  Results  and  Analysis 


4-1  Overview 

This  chapter  presents  performance  results  of  a  symbol  timing  dither  system  un¬ 
der  multipath  conditions.  Section  4.2  presents  a  performance  comparison  of 
differing  (3  valued  SRRC  pulse  shapes  to  determine  what  effect  the  severity  of  the 
pulse  shape  roll  off  has  on  bit  error  performance.  Section  4.4  presents  results  of  a  two 
ray  multipath  model  under  a  variety  of  symbol  dither  conditions,  and  Section  4.5  ex¬ 
amines  the  performance  of  a  four  ray  multipath  model.  Both  models  are  constructed 
according  to  the  parameters  provided  in  Section  3.4.1. 

Results  in  this  chapter  were  generated  using  1V=1000  randomly  generated  sym¬ 
bols  for  both  training  of  the  non-cooperative  CMA  receiver  and  a  unique  random  set 
of  1000  symbols  for  symbol  estimation,  with  N  symbols  being  transmitted  over  each 
of  the  random  channels.  The  BER  results  are  then  averaged  together  to  provide  a 
robust  picture  of  system  performance  under  the  tested  conditions.  The  confidence 
intervals  shown  in  figures  within  Sections  4.4  and  4.5  correspond  to  95%  confidence 
intervals  based  on  (4.1)  derived  from  [7] 


3hZbL 


a 


n 


(4.1) 


where  Z°l  corresponds  to  1.96  for  a  95%  confidence,  n  is  the  number  of  samples  (in 
this  case  the  number  of  randomly  generated  channels),  and  cr  is  the  standard  deviation 
of  the  calculated  Pb  data.  The  standard  deviation  is  calculated  by  taking  the  BER 
of  each  channel  and  then  taking  the  standard  deviation  of  the  resulting  vector.  The 
results  presented  in  Section  4.4  and  Section  4.5  are  separated  between  well-behaved 
channel  performance  plots  and  poorly  performing  channel  plots.  The  process  used  to 
define  a  channel  as  well-behaved  or  poorly  performing  is  presented  in  Section  4.3,  and 
results  in  35  well-behaved  channels  and  15  poor  performance  channels  for  the  two  ray 
multipath  model  and  44  well-behaved,  6  poorly  performing  channels  for  the  four  ray 
multipath  model.  Table  4.1  provides  an  overview  of  the  results  presented. 
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Table  4.1:  Table  of  experiments 


Figure 

P 

Max  Dither  ( Tsyrn ) 

Eb/N0 

Multipath  Model 

4.1 

0 

full  range 

9  dB 

none 

4.2 

0.25 

full  range 

9  dB 

none 

4.3 

0.5 

full  range 

9  dB 

none 

4.8 

0.5 

0.6 

3-15  dB 

random  2  ray 

4.9 

0.5 

0.99 

3-15  clB 

random  2  ray 

4.10 

0.5 

full  range 

9  dB 

random  2  ray 

4.11 

0.5 

0.6 

3-9  dB 

bad  random  2  ray 

4.12 

0.5 

0.99 

3-15  dB 

random  4  ray 

4.13 

0.5 

0.6 

3-15  dB 

random  4  ray 

4.14 

0.5 

full  range 

9  dB 

random  4  ray 

4.15 

0.5 

0 

3-9  dB 

bad  random  4  ray 

4.16 

0.5 

0.6 

3-9  dB 

bad  random  4  ray 

4.17 

0.5 

0 

9  dB 

bad  random  4  ray 

4-  2  Pulse  Shape  Performance 

This  section  examines  the  effects  of  varying  the  roll  off  factor  (3  of  the  SRRC 
pulse  shape  used  throughout  this  thesis.  Figure  4.1,  Figure  4.2,  and  Figure  4.3  show 
the  Pb  achieved  by  each  equalizer  structure  as  the  timing  symbol  dither  range  is  varied 
in  the  absence  of  any  multipath  interference.  All  plots  shown  are  held  constant  at 
Eb/Na  of  9  dB  for  comparison.  As  the  (3  value  increases  from  0  to  0.5,  the  pulse  shape 
contains  less  energy  outside  the  main  lobe  of  the  pulse,  resulting  in  less  ISI  as  the 
symbol  dither  increases.  Also,  for  (3= 0.5,  the  cooperative  MSE  equalizer  structure 
suffers  no  significant  degradation  in  terms  of  bit  error  performance  until  the  symbol 
dither  range  exceeds  0.6Tsym.  This  is  a  marked  increase  from  the  (3=0  case,  where 
the  MSE  equalizer  begins  to  degrade  at  a  symbol  dither  range  of  0.2 Tsym.  The  results 
for  (3= 0.25  fall  somewhere  between  the  performance  of  (3=0  and  (3= 0.5,  offering  an 
intermediate  step  that  generates  less  ISI  than  the  (3=0  case. 
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Figure  4.1:  Bit  Error  Comparison  for  SRRC  j3— 0  at  Eb/N0  =  9  dB 


Figure  4.2:  Bit  Error  Comparison  for  SRRC  /3=0.25  at  Eb/N0  =  9  dB 


For  no  symbol  timing  dither,  all  three  pulse  shapes  and  equalizers  result  in  Pb 
=  10“3  or  better,  and  are  close  to  the  analytical  curve,  as  is  expected  from  validation 
exercises  shown  in  Chapter  III.  As  /3  increases,  there  is  a  small  improvement  across 


all  equalizer  structures  as  the  symbol  timing  dither  range  increases;  10-20%  of  an 
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Figure  4.3:  Bit  Error  Comparison  for  SRRC  /?=0.5  at  Eb/Na  =  9  dB 

order  of  magnitude.  As  f3= 0.5  gives  slightly  better  bit  error  performance  across  all 
equalizer  structures,  all  other  results  presented  in  this  chapter  utilize  this  pulse  shape. 

4-3  Multipath  Channel  Validation 

To  validate  each  randomly  generated  multipath  channel,  the  BER  was  recorded 
over  5000  transmitted  symbols.  These  bit  error  rates  were  then  used  to  sort  the 
channels  from  best  to  worst  in  terms  of  Pb  performance.  For  both  the  two  ray  and 
four  ray  groups  of  multipath  channels,  an  order  of  magnitude  breakpoint  existed 
in  channel  performance.  Channel  models  that  provided  equalization  results  better 
than  this  threshold  are  defined  as  well-behaved,  and  those  that  generated  BERs  in 
excess  of  this  threshold  are  defined  as  poorly  performing.  Examples  of  these  randomly 
generated  multipath  channels  are  shown  in  Figure  4.4  and  Figure  4.5. 

Figure  4.6  shows  the  results  of  applying  the  determined  MSE  equalization  filter 
to  a  well-behaved  four  ray  multipath  channel  with  no  timing  symbol  dither.  The 
large  centered  impulse  demonstrates  that  a  majority  of  the  multipath  effects  have 
been  removed  by  the  equalization  process.  Figure  4.7  shows  the  results  of  this  process 
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Figure  4.4:  Poor  Performance  Two  Ray  Multipath  Model  Normalized  channel  be¬ 
fore  equalization  at  Eb/N0  =  9  dB 


Figure  4.5:  Good  Performance  Four  Ray  Multipath  Model  Normalized  channel 

before  equalization  at  Eb/N0  =  9  dB 


for  a  poorly  performing  two  ray  multipath  model  with  no  timing  symbol  dither.  The 
lack  of  a  large  impulse  relative  to  the  others  indicates  poor  equalization  performance. 
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Figure  4.6:  Good  Performance  Four  Ray  Multipath  Model  Normalized  channel  after 
equalization  at  Eb/N0  =  9  dB 


Figure  4.7:  Poor  Performance  Two  Ray  Multipath  Model  Normalized  channel  after 
equalization  at  Eb/N0  =  9  dB 

Inspection  of  the  coefficients  of  the  poorly  performing  channels  in  terms  of  Pb  for 
both  the  2  ray  and  4  ray  models  revealed  that  each  contained  a  multipath  reflection 
of  near  equal  normalized  power  as  the  line  of  sight  path.  In  order  to  understand 
this  behavior,  the  z  transform  of  a  poorly  performing  channel  is  examined.  The  z 
transform  is  given  in  (4.2)  [11]. 

OO 

x(z)  =  J2xin}z~n  (4-2) 
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The  two  ray  channel  model  for  a  poor  performance  channel  can  be  approximated  as 
h  =  [1  1],  Taking  the  z  transform  of  h  yields  the  result,  shown  in  (4.3). 

H(z)  =  l  +  z~1  =  (::  +  1)  (4.3) 

z 

This  results  in  a  zero  on  the  unit  circle,  and  as  the  equalizer  attempts  to  invert 
the  channel  effects  to  remove  them,  this  zero  will  become  a  pole  after  inversion.  This 
pole  on  the  unit  circle  means  that  the  channel  approximated  by  h  above  does  not  invert 
well,  and  prevents  proper  calculation  of  equalization  filter  coefficients.  These  poorly 
inverting  channels  dominate  the  Pb  performance  when  all  50  random  channels  are 
averaged  together,  so  the  poor  performance  channels  are  separated  out  and  presented 
separately  from  the  well-behaved  channel  performance  when  presented  in  Chapter  IV. 

4-4  Two  Ray  Model  Multipath  Performance 

The  results  presented  in  this  section  are  divided  into  two  parts,  with  the  perfor¬ 
mance  over  well-behaved  two  ray  random  multipath  models  presented  first,  followed 
by  results  achieved  using  poorly  behaved  two  ray  models  presented  in  Section  4.4.1. 
Figure  4.8  and  Figure  4.9  present  the  results  achieved  using  the  randomly  generated 
well-behaved  two  ray  multipath  channels.  For  a  symbol  timing  dither  range  of  0.6 Tsym, 
the  MSE  cooperative  receiver  achieves  a  BER  of  Pb  =  2  x  10~2  with  approximately 
7  dB  less  Eb/N0.  Through  Eb/N0  =  15  dB,  both  non-cooperative  receiver  structures 
BERs  remain  in  excess  of  1 0  ~  2 .  The  BER  performance  of  the  MSE  cooperative  re¬ 
ceiver  asymptotically  approaches  Pb  =  5  x  1CT3  at  Eb/N0  =  11  dB,  and  remains  stable 
for  higher  values  of  Eb/N0.  Due  to  the  power  in  the  multipath  model  being  normal¬ 
ized,  increasing  the  signal  power  also  increases  the  power  of  the  multipath  reflection, 
creating  a  BER  breakpoint  for  the  MSE  equalizer  structure  used  in  these  simulations. 
Figure  4.9  also  shows  this  asymptotic  behavior  in  the  MSE  equalizer  performance, 
but  the  BER  achieved  at  Eb/Na  has  yet  to  surpass  the  breakpoint  value  of  5  x  10”3 
noted  in  Figure  4.8.  For  the  full  range  of  symbol  dither  presented  here,  neither  of 
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Figure  4.8:  Two  Ray  Random  Multipath,  Good  Channels,  0.6 Tsym  Max  Delay 

the  non-cooperative  equalizers  achieve  BERs  lower  than  approximately  0.1,  while  the 
MSE  equalizer  achieves  p,  values  an  order  of  magnitude  better  (p,  =  10-2).  Fig¬ 
ure  4.10  shows  statistically  significant  performance  improvement  of  the  cooperative 
MSE  equalizer  over  the  non-cooperative  equalizers  explored  in  this  thesis  over  the  full 
range  of  possible  symbol  timing  dither  values. 
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Figure  4.9:  Two  Ray  Random  Multipath,  Good  Channels,  0.99 Tsym  Max  Delay 

4-4-1  Poorly  Performing  Two  Ray  Model  Performance.  Figure  4.11  shows 
the  performance  achieved  using  the  15  poorly  behaved  random  two  ray  multipath 
channels.  The  difficulty  in  inverting  the  channel  for  equalization  discussed  in  Sec¬ 
tion  3.4.1  results  in  the  multipath  dominating  bit  error  performance,  with  all  equal¬ 
izer  structures  generating  BERs  in  excess  of  P, J,  =  10"1  for  a  moderate  symbol  timing 
dither  range  value  of  0.6 Tsym.  This  behavior  indicates  that  a  system  using  the  dither 
algorithm  method  equalizer  techniques  presented  is  sensitive  to  the  multipath  condi¬ 
tions  present  in  the  channel  at  the  time  of  transmission.  This  is  not  the  fault  of  the 
dither  algorithm,  however,  as  Figure  4.10  shows  a  near  order  of  magnitude  perfor¬ 
mance  degradation  of  the  CMA  equalizer  as  compared  to  the  MMSE  equalizer  even 
for  no  symbol  dither  conditions. 
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Figure  4.10:  Two  Ray  Random  Multipath,  Good  Channels,  Eb/N0  =  9  dB 


Eb/N0  (dB) 


Figure  4.11:  Two  Ray  Random  Multipath,  Bad  Channels,  0.6 Tsym  Max  Delay 


4-5  Four  Ray  Model  Multipath  Performance 


Figure  4.12  and  Figure  4.13  characterize  the  BER  behavior  of  the  symbol  timing 


dither  system  in  the  presence  of  the  44  random  four  ray  multipath  channels.  The  same 
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asymptotic  behavior  noted  for  the  two  ray  multipath  models  occurs  for  the  four  ray 
multipath  as  well,  occurring  at  a  BER  of  approximately  Pb  =  1  x  10"4  for  Eb/N0  = 
13  dB  for  the  MSE  equalizer  as  seen  in  Figure  4.13.  The  MSE  equalizer  provides  a 
BER  performance  gain  of  approximately  7  dB  relative  to  explored  non-cooperative 
methods  at  Pb  =  2x  10~2  and  0.6 Tsym  symbol  timing  dither  range.  Using  the  maximum 
allowable  symbol  timing  dither  of  0.99 Tsym  (Figure  4.12),  the  MSE  equalizer  provides 
6  dB  of  gain  over  all  non-cooperative  methods  at  a  BER  of  10_1.  At  Eb/N0  =  15  dB, 
the  MSE  equalizer  outperforms  the  non-cooperative  equalizers  by  nearly  two  orders 
of  magnitude,  with  BERs  of  2  x  10-3.  Also,  the  non-cooperative  equalizers  produce  a 
BER  in  excess  of  10_1  for  all  Eb/N0  values  examined,  indicating  a  theoretical  limit  for 
non-cooperative  performance  under  the  tested  multipath  and  symbol  timing  dither 
conditions.  Figure  4.14  shows  an  approximate  order  of  magnitude  BER  performance 
improvement  using  a  cooperative  MSE  equalizer  attuned  to  symbol  timing  dither  over 
non-cooperative  equalization  methods  for  all  possible  ranges  of  symbol  timing  dither. 


Figure  4.12:  Four  Ray  Random  Multipath,  Good  Channels,  0.99 Tsym  Max  Delay 
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Figure  4.13:  Four  Ray  Random  Multipath,  Good  Channels,  0.6 Tsym  Max  Delay 

4-5.1  Poorly  Performing  Four  Ray  Model  Performance.  Figure  4.15  and 
Figure  4.16  show  the  performance  achieved  using  the  6  poorly  behaved  random  four 
ray  multipath  channels.  The  difficulty  in  inverting  the  channel  for  equalization  dis¬ 
cussed  in  Section  3.4.1  results  in  the  multipath  dominating  bit  error  performance  as 
seen  for  the  poorly  performing  two  ray  model  channels,  with  all  equalizer  structures 
generating  BERs  in  excess  of  P, J,  =  ICE1  for  a  both  a  symbol  timing  dither  range  value 
of  0.4  and  0.6 Tsym.  This  behavior  indicates  that  the  dither  algorithm  method  utiliz¬ 
ing  the  equalizer  techniques  presented  is  sensitive  to  the  multipath  conditions  present 
in  the  channel  at  the  time  of  transmission.  As  in  the  two  ray  model  results,  this  is 
not  the  fault  of  the  dither  algorithm.  Figure  4.14  shows  a  near  order  of  magnitude 
performance  degradation  of  the  CMA  equalizer  as  compared  to  the  MMSE  equalizer 
for  no  symbol  dither  conditions.  Figure  4.17  repeats  the  experiments  performed  in 
Figure  4.15  and  Figure  4.16  with  no  symbol  timing  dither,  resulting  in  comparable 
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Figure  4.14:  Four  Ray  Random  Multipath,  Good  Channels,  Eb/N0= 9  dB 


bit  error  performance  for  all  equalization  schemes.  This  indicates  that  the  sensitivity 
of  the  symbol  timing  dither  system  to  multipath  conditions  noted  above  is  due  to  the 
difficulty  involved  in  equalizing  certain  multipath  channels,  not  any  internal  flaw  with 
the  dither  scheme. 
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CL 


Figure  4.15: 


JD 

CL 


Figure  4.16: 


Eb/Na  (dB) 


Four  Ray  Random  Multipath,  Bad  Channels,  0ATsym  Max  Delay 


Eb/N0  (dB) 


Four  Ray  Random  Multipath,  Bad  Channels,  0.6 Tsym  Max  Delay 
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CL 


Figure  4.17: 


Eb/N0  (dB) 


Four  Ray  Random  Multipath,  Bad  Channels,  no  symbol  timing  dither 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions 

This  thesis  tests  the  performance  of  an  LPI  system  implemented  using  the  sym¬ 
bol  timing  dither  in  the  presence  of  multipath  channel  interference  [8].  A  coopera¬ 
tive  receiver  structure  is  shown  to  provide  significant  performance  improvement  over 
tested  non-cooperative  receiver  structures  for  a  binary  DPSK  modulated  waveform. 
The  timing  dither  coupled  with  multipath  interference  creates  ISI  that  is  not  explicitly 
accounted  for  in  non-cooperative  receiver  designs,  with  the  maximum  allowed  symbol 
timing  dither  range  and  pulse  shape  used  determining  the  severity  of  the  induced  ISI 
for  a  particular  multipath  interference  channel. 

For  both  the  two  and  four  ray  randomly  generated  multipath  channel  models, 
a  cooperative  MMSE  equalizer  structure  having  knowledge  of  and  synchronization 
with  the  transmitted  symbol  timing  dither  can  be  used  to  achieve  a  factor  of  10  BER 
performance  gain  over  non-cooperative  equalizer  methods.  Typical  gains  of  the  MSE 
equalizer  over  tested  non-cooperative  equalization  techniques  of  6-8  dB  were  achieved 
at  BERs  of  1),  =  10-1.  When  using  the  maximum  allowable  range  of  symbol  tim¬ 
ing  dither,  0.99 Tsym,  non-cooperative  equalization  techniques  consistently  produced 
BERs  in  excess  of  10_1  even  for  large  values  of  Eb/Na.  The  MSE  equalizer  per¬ 
formance  also  exhibited  an  upper  limit,  asymptotically  approaching  BERs  of  P&  = 
5  x  10“3  and  1  x  10-4  for  the  two  and  four  ray  multipath  models,  respectively.  The 
MMSE  equalizer  structure  used  in  this  thesis  is  computationally  expensive,  requiring 
a  matrix  inverse  operation  for  each  transmitted  symbol.  Also,  for  certain  randomly 
generated  multipath  channels,  this  inversion  operation  was  not  well-behaved,  and  re¬ 
sulted  in  poor  equalization  results  for  all  equalizers.  This  is  not  necessarily  due  to 
equalizer  or  system  issues,  as  some  channels  are  simply  difficult  to  equalize  effectively. 
This  indicates  a  certain  amount  of  sensitivity  to  the  multipath  channel  conditions  on 
the  part  of  a  symbol  timing  dither  system.  Under  well-behaved  and  invertible  multi- 
path  interference  conditions,  effective  symbol  recovery  of  a  transmitted  LPI  waveform 
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making  use  of  symbol  timing  dither  is  feasible  assuming  a  cooperative  receiver  with 
knowledge  of  the  dither  pattern. 

5.2  Recommendations  for  Further  Study 

Several  areas  for  further  study  related  to  this  thesis  are  summarized  in  this 
section. 

It  is  possible  to  alter  the  design  of  the  symbol  timing  dither  algorithm.  The 
algorithm  used  for  this  thesis  always  followed  a  uniform  probability  density  function, 
but  many  other  distributions  could  be  used.  Also,  the  results  here  are  limited  to 
binary  DPSK  modulation.  Other  modulation  techniques  could  be  investigated  for  use 
with  an  LPI  system  based  on  symbol  timing  dither. 

A  limited  set  of  possible  multipath  models  were  used  in  this  thesis  to  characterize 
a  symbol  timing  dither  waveform’s  performance  in  multipath  interference.  These 
models  used  fixed  time  delays  to  ensure  that  all  the  channels  used  were  linear  time 
invariant.  Further  research  could  characterize  the  symbol  timing  dither  waveform’s 
performance  under  time  varying  multipath  conditions. 

Additional  cooperative  and  non-cooperative  receiver  designs  could  be  examined 
for  use  with  symbol  timing  dither.  A  small  subset  of  the  possible  receiver  designs  were 
tested  in  this  thesis,  and  these  could  be  improved  upon  or  new  receiver  structures  and 
equalization  techniques  could  be  examined. 
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